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ABSTRACT 

This  project  deals  with  two  types  of  nonlinear  inter¬ 
action  of  intense  laser  beam  with  dielectric  medium.  The  first 
and  extensively  studied  phenomenon  is  the  nonlinear  dc  polari¬ 
zation  that  is  caused  by  the  second  order  nonlinear  polarization. 
The  quartz  crystal  is  chosen  as  the  dielectric  medium.  A 
quantitative  relationship  between  the  dc  polarization  and  the 
intensity  of  the  propagating  laser  beam  is  given.  A  convenient 
method  of  detecting  dc  polarization  is  presented.  With  this 
technique  the  elements  of  the  second  order  polarization  co¬ 
efficient  can  be  determined  experimentally.  By  considering  a 
suitable  detecting  system  wioh  a  convenient  configuration  of 
the  quartz  crystal,  it  is  shown  that  the  output  voltage  of  the 
detector  is  linearly  proportional  to  the  intensity  of  the  laser 
pulse.  Thus,  the  possibility  of  using  this  principle  to  build 
a  transmission  type  of  meter  for  measuring  power  in  high-power 
laser  pulses  is  presented.  Preliminary  experimental  results  on 
the  correlation  between  the  spikings  in  the  dc  polarization  and 
in  the  fundamental  frequency  of  laser  oscillation  are  presented. 

The  second  type  of  nonlinear  polarization  presented 
in  this  report  is  the  initiation  of  a  study  on  the  generation  of 
subharmonics  at  optical  frequencies  using  high  intensity  laser 
sources.  The  physical  mechanism  that  is  utilized  for  this 


purpose  is  the  nonlinear  interaction  of  electromagnetic 


iii 


radiation  with  certain  types  of  crystalline  media  such  as  quartz 
and  KDP.  This  nonlinear  interaction  results  in  parametric  sub¬ 
harmonic  generation  under  favorable  conditions.  The  experimental 
arrangement  used  for  this  investigation  is  outlined. 


1.  INTRODUCTION  I  f 


The  purpose  of  this  project  was  to  study  the  nonlinear 
polarization  effects  in  crystals  when  high  intensity  electro¬ 
magnetic  wa,-e  propagates  through  them.  This  was  made  possible 
by  the  availability  of  high  power  lasers.  Two  types  of  non- 
linearities  are  treated  in  this  report.  They  are: 

1.  Second  order  nonlinear  d.c.  polarization. 

2.  Generation  of  subharmonics  due  to  nonlinear  polarization. 

The  earliest  result  of  the  efforts  in  the  investigation 
of  nonlinear  polarization  using  laser  beams  was  reported  by 
Franken,  et .  al.1.  They  observed  second  harmonic  generation 
by  passing  ruby  laser  beam  through  crystalline  quartz  medium. 

The  production  of  second  harmonic  can  be  explained  by  considering 
the  following  simple  scalar  mathematical  model  relating  the 
polarization  p  and  the  electric  field  intensity  E. 

P  =  ajE  +  a2E2  +  a^3  +  .  (l) 

Here  a^,  a^,  a^,  ........  are  called  the  first,  second,  third, 

.  order  polarization  coefficients.  In  writing  Eq.  (l) 

I 

the  gradients  of  the  E  field  have  been  neglected  for  simplicity. 

If  the  electric  field  is  caused  by  the  laser  beam,  E  can  be 

written  as  E  =  E  cos  oot,  cc  being  the  angular  frequency  of  the 
1  o 

I 

radiation  field.  Substituting  this  value  of  E  in  Eq.  (l),  one 
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p  =  a,E„  con  o)t  +  — -  ( 1  +  cos  2o.'t) 

io  2  o 


a  n 

— ^  (3  cos  cut  +  cos  3cot)  + 


(2) 


It  can  be  readily  observed  that  the  secc  d  term  in  the  right  hand 

side  of  Eq.  (2)  generates  the  second  harmonic.  In  general,  the 

coefficients  a^,  a^,  a^,  ......  are  tensors.  However,  a  dis- 

p 

cussion  of  their  tensor  form  is  deferred.  Franken  and  Ward 
have  shown  in  detail  that  the  coefficient  a^  will  be  nonzero 
only  in  those  media  which  posses  lack  of  inversion  Symmetry. 

From  Eq.  (2)  it  can  be  seen  that  the  second  harmonic 
generation  is  accompanied  by  a  d.c.  component.  Thus,  the 
successful  observation  of  the  second  harmonic  output  initiated 

3 

interest  in  investigating  the  d.c.  polarization.  Bass,  et.  al. 

have  reported  observing  the  d.c.  polarization  in  KDP  and  KD^P 

4 

crystals.  The  authors  have  applied  this  phenomenon  in  quartz 
crystal  for  laser  power  measurement.  Section  2  of  this  report 
presents  the  details  of  the  extensive  work  that  was  performed  on 
the  d.c.  polarization. 

The  preliminary  investigation  made  to  generate  sub¬ 
harmonics  using  high  intensity  laser  sources  is  presented  in 

Section  3*  The  subject  of  subharmonic  oscillation  has  been 

2  5  6  5 

dealt  with  by  various  authors .  Armstrong,  et.al.  ,  and 

2 

Franken  and  Ward  have  considered  the  propagation  of  three 
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signals  through  a  nonlinear  medium.  Under  degenerate  conditions, 
energy  is  transferred  f^om  a  signal  at  frequency  2<d  to  the  signal 
at  frequency  (o.  Kingston  and  McWhorter  have  analyzed  the  problem 
using  coupled  mode  theory.  Up  to  the  present  time,  to  the  best  of 
the  authors*  knowledge,  no  experimental  evidence  of  the  existence 
of  optical  subharmonics  has  been  demonstrated,  and  our  efforts  are 
aimed  toward  this  goal. 


4 

2.  NONLINEAR  D.C.  POLAI  ..ZATION  AT  OPTICAL  FREQUENCIES 


* 


2.1  Introduction 

The  generalized  expression  for  polarization,  in  the 
absence  of  any  external  biasing  electric  or  magnetic  field,  may 
be  written  as  a  power  series  in  terms  of  the  components  of  the 
electric  field  intensity  and  its  gradient. 


>,  =  e  : 

1  0  L 


xuEj +  Ek +  xljk  EA  +  xij^  Ej  uk  Ei 


+  xi^l  Ej  \  E« + 


(3) 


J. 

where  pi  is  the  i  spatial  component  of  polarization  p,  X*s 
represent  the  various  orders  of  polarization  coefficient  tensors, 
V's  denote  the  gradient  operation,  E‘s  are  the  spatial  components 
of  the  electric  field  intensity  and  eQ  the  free  space  permlltivity . 

O 

Franken  and  Ward  have  discussed  the  physical  significance  of  the 
various  terms  in  Eq.  (3).  It  is  enough  to  mention  here  that  the 
only  term  in  the  right  hand  side  that  can  cause  d.c.  polarization 
in  the  medium  is  the  second  order  term  E^  E^.  The  magnitude 

of  d.c.  polarization  due  to  terms  of  higher  order  than  those  in 
Eq*  (3)  can  be  neglected  as  compared  with  that  of  the  second  order 
term  since  the  contribution  to  polarization  decreases  as  the  ratio 
of  the  electric  field  intensity  E  of  the  electromagnetic  wave  to 
the  atomic  electric  field  intensity  Ea^omj_c  (E/^atomic^  ~or  each 

p 

additional  E  factor  added  . 
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Retaining  only  the  second  order  term,  Eq.  (3)  can  be 


written  as 


p.  =  X.  n  E  .  E, 
i  ijk  j  1: 


where  ^  is  a  third  rank  tensor.  It  will  oe  shown  in  Section  2.9 

that  the  second  order  polarization  tensor  is  related  to  the  linear 

a 

electro-optic  tensor.  The  experimental  results  of  Bass,  et.  al. 
on  KDP  and  KD^P  crystals  and  one  of  the  results  presented  in  this 
report  confirms  the  above  theory. 


This  section  deals  with  the  observation  of  the  d.c. 


polarization  in  quartz  crystal  and  by  measurement  of  this  d.c. 

polarization,  calculating  one  of  the  second  order  nonlinear 

polarization  coefficients.  The  theory  of  the  propagation  of 

the  laser  beam  along  the  optic  axis  of  the  quartz  dielectric 

4 

medium  has  been  treated  elsewhere  and  a  brief  review  is  given 
in  Section  2.2.  The  d.c.  polarization  developed  by  the  laser 
beam  is  detected  due  to  its  interaction  with  an  external  detecting 
system.  This  is  explained  in  Section  2.3  by  assuming  a  simplified 
model.  The  equivalent  circuit  representation  is  then  derived  for 
the  system.  Prom  both  the  circuitry  and  field  theory  aspects,  it 
is  shown  that  no  energy  conversion  is  possible  from  optical 


frequency  to  d.c.  by  making  use  of  this  phenomenon.  In  Sections 
2.4  and  2.3,  the  theoretical  and  constructional  details  of  the 
quartz  detector  are  presented.  Tb  experimental  arrangement  and 
results  on  quartz  crystal  are  given  in  Section  2.6.  The  results 
of  Sections  2.7,  2.8  and  2.9  include  the  observation  of  d.c. 
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polarization,  its  angular  dependence  with  respect  to  crystal  axes 
orientation  and  the  estimation  of  the  second  order  nonlinear  polar¬ 
ization  coefficient.  Section  2.10  gives  the  experimental  results 
of  the  spiking  phenomenon  in  the  d.c.  polarization.  The  results 
are  summarized  in  Section  2.11. 

2 . 2  Propagation  in  a_  Infinite  Medium 

The  propagation  of  an  electromagnetic  wave  through  a 
nonlinear  medium  has  already  been  discussed  in  detail  by  various 
authors-^  from  both  the  quantum  mechanical  and  phenomenological 
approaches.  Vrhen  considering  the  interaction  between  various  waves 
propagating  in  a  medium,  one  has  to  take  into  account  such  effects 
as  the  dispersion  in  the  medium,  the  phase  matching  of  the  various 
components  and  so  forth.  However,  in  the  case  of  d.c.,  the  phase 
velocity  is  infinite  and  there  is  no  propagation  of  wave  at  zero 
frequency.  An  analysis  for  the  d.c.  case  has  been  presented  by 
the  authors  elsewhere  and  only  a  brief  summary  is  presented  here. 

The  nonlinear  dielectric  considered  is  crystalline 
quartz.  The  medium  is  assumed  to  be  infinite  in  extent  and  non- 
dissipative.  The  second  order  polarization  tensor  and  the 
piezoelectric  tensor,  both  of  rank  3>  are  determined  by  the 
same  symmetry  conditions  and  hence  have  the  same  form.  Thus, 

Q 

the  following  contracted  form  of  the  piezoelectric  notation 
can  be  used  to  write  in  Eq.  (4). 


"i  jk  Aij 


(i  =  1,  2,  3;  J  =  1,  2,  .  ,6)  (3) 


. . . . 
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For  quartz  crystal  belonging  to  class  32,  xu  can  be  expressed  In 
the  form  similar  to  its  piezoelectric  tensor^  which  is  shown  below 

/a  -a  0  0  0  0  \ 


XU  e0 


0  0  0 
0  -0  -2a 


From  Eqs.  (4),  (5),  and  (6) 


»x  -  °<Ex2  -  Ey2)  +  eEyEz 


Py  -  '  @VZ  '  31  ExEy 


n_  =  0 


The  special  case  of  propagatJon  along  the  optic  axis 

will  now  be  considered.  The  orientation  of  the  crystal  axes  and 

the  direction  of  polarization  of  the  beam  are  shown  in  Figure  1. 

The  laser  beam  is  assumed  to  be  linearly  polarized  in  a  direction 

making  an  angle  8  with  the  x-axis  of  the  crystal.  In  Figure  1, 

x*x*  and  yty*  are  arbitrary  choices  of  coordinate  axes,  and  xx 

.4 

and  yy  are  the  cyrstal  axes.  Then,  it  can  be  shown  that  the 
d.c.  polarization  along  the  two  crystal  axes  are  given  by 


—  V 
Px  2 


cos  26 


Py  "  2 


sin  26 


where  E  =  Z  cos  o>t  is  the  electric  field  intensity  of  the  incident 
radiation.  Figure  1  along  with  Eqs  (8)  and  (9)  describe  the 
angular  dependence  of  d.c.  polarization.  It  can  be  observed  that 
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as  the  direction  of  polarization  of  the  incident  radiation  is 
rotated  through  an  angle  6,  the  d.c.  polarization  is  rotated 
through  an  angle  20  in  the  opposite  direction. 

2 . 3  Detecting  Technique  and  Circuit  Considerations 

This  section  is  devoted  to  the  technique  by  which  the 
d.c.  polarization  is  detected  and  an  equivalent  circuit  model  is 
developed  for  the  detecting  system.  In  order  to  develop  the  con¬ 
cept,  a  simple  model  is  chosen  and  an  analysis  of  it  is  made. 

The  nonlinear  polarization  can  be  detected  with  the  use 
of  an  external  circuit  arrangement  of  the  type  described  below. 
The  medium  in  which  the  nonlinear  polarization  is  established  is 
made  part  of  a  capacitor  that  is  formed  by  two  electrodes  that 
are  placed  on  opposite  sides  of  the  dielectric.  The  capacitor 
thus  formed  is  then  connected  to  an  external  detector.  The 
arrangement  shown  In  Figure  2  represents  the  simple  case  of  a 
parallel  plate  capacitor  formed  with  the  nonlinear  dielectric 
medium. 


The  voltage  Vq  across  the  capacitor,  neglecting  the 
fringe-effects  is  given,  by 


2P0d 

e 


(10) 


where  PQ  is  the  uniform  polarization  induced  by  the  laser  beam 
in  the  entire  dielectric  medium  along  y-direction,  2d  the  thickness 
of  the  dielectric  and  is  also  the  distance  between  capacitor  plates, 
and  e  the  permittivity  of  the  dielectric  at  low  frequencies. 
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The  capacitor  shown  in  Figure  2  is  now  connected  to  an 
external  detector  with  an  input  capacitance  C  and  an  input  re¬ 
sistance  R.  The  nodal  equation  at  node  A  is 

^  dV 

■a§  +  ~R  +  c~dt  =  0  (11) 

where  VQ  is  the  instantaneous  voltage  across  the  capacitor  plates 
and  +Q  and  -Q  are  the  charges  on  them.  The  charge  Q  on  each  plate 
is  due  to  the  external  polarization  source  plus  the  depolarizing 
effect  in  the  material.  Assuming  unit  area  for  the  capacitor  plate, 
the  total  charge  on  the  capacitor  is 


Q  =  £E  +  Pq  (12) 

where  E  is  the  electric  field  intensity  in  the  capacitor.  Sub¬ 
stituting  Eq.  (12)  in  Eq.  (11),  one  has 


dVo 


(13) 


where  E  has  been  replaced  by  V^/^d.  For  the  nonfringing  case  that 
is  under  consideration,  e/2d  is  the  capacitance  of  the  capacitor 
formed  with  the  nonlinear  dielectric.  If  this  capacitance  is 
designated  by  Cq  then  Eq.  (13)  can  be  rewritten  as 


dV. 


V. 


dP. 


dt  +  (CQ+C)  R  +  (CQ+C)  '  dt 


=  0 


Now,  PQ  can  be  related  to  Vq  from  Eq.  (10)  by 


Po  -  V 


Q  C4 


(14) 


(15) 


10 


where  Is  an  equivalent  capacitance  and  is  given  by 

C4  =  2d  "  CQ 


(16) 


Although  and  are  equal  in  the  present  case,  it  need  not  be  the 
case  always.  If,  for  example,  there  is  an  air  gap  between  the  plate 
and  the  dielectric,  and  Cq  will  be  different.  One  has  to  dis¬ 
tinguish  this  difference  in  an  actual  system.  From  Eqs.  (l4)  and 
(15),  the  equivalent  circuit  for  the  system  shown  in  Figure  3  can 
be  synthesized  and  is  given  in  Figure  4.  It  can  be  seen  that  the 

available  polarization  source  voltage  V  is  less  than  Vn,  since 

e 

<  Cq  and  is  given  by 


V  -  y  — 
e  v«  c« 


(17) 


The  output  response  to  a  continuous  laser  beam  that  is 
turned  on  at  time  t  =  0,  is  given  in  Figure  5.  It  is  seen  that  the 
output  voltage  reduces  to  zero  exponentially.  Thus,  no  continuous 
d.c.  power  is  delivered  to  the  output  except  during  the  transient 
condition;  that  is,  no  optical  power  rectification  is  possible. 

This  could  al3o  be  explained  from  the  field  theory  point  of  view. 
The  network  done  by  the  field  on  the  system  is  given  by  flT  •  -^3. 
For  d.c.  polarization  t—  =  0  and  hence  there  is  nr  net  energy 
transfer  in  the  steady  state  condition. 

2.4  Theory  on  the  Quartz  Detector 

In  this  section  a  suitable  configuration  of  electrodes 
is  determined  for  the  type  of  the  detector  described  in  the 
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previous  section.  Consider  a  circular  cylindrical  quartz  rod 
with  the  optic  axis  oriented  along  the  length  of  the  rod.  Let 
the  laser  beam  propagate  along  the  axial  direction.  Figure  6 
shows  the  cross-section  of  the  cylindrical  quartz  rod  of  radius 
b  with  Its  optic  axis  perpendicular  to  the  plane  of  the  paper.  The 
incident  beam  is  assumed  to  be  linearly  polarized  and  cylindrical 
in  cross-section.  The  radius  of  the  beam  is  a. 

Let  the  laser  beam  cause  a  uniform  d.c.  polarization 

given  by  Eqs.  (8)  and  (9)  in  the  transverse  direction  making  an 

angle  0  with  the  axis  of  the  crystal.  Since  quartz  is  a  uniaxial 

crystal,  the  x-y  plane  is  isotropic.  This  can  be  considered  as  a 

two-dimensional  electrostatic  problem  which  has  already  been 

4 

solved  .  The  potential  outside  the  crystal  is  given  by 

V  -  ,  e  -a?e-  ■+  Tf  7  «M)  (18) 

where  a  is  the  nonlinear  polarization  coefficient  described  in 
Eq.  (6),  TJ  is  the  intrinsic  impedance  of  the  quartz  medium,  eQ  is 
the  free-space  permittivity,  er  the  relative  dielectric  constant 
of  quartz  at  low  frequencies  and  is  the  in*ensity  of  the  laser 
beam.  The  equipotential  lines  described  by  Eq.  (l8)  with  0  =  0, 
are  shown  in  Figure  7.  It  can  be  observed  that  the  equation  r  =  k 
cos  0.  where  k  is  an  arbitrary  constant,  describes  equipotential 
surfaces  which  are  pairs  of  circles.  For  the  construction  of  the 
detector,  one  such  pair  corresponding  to  r  =  k^  cos  0  is  chosen. 
Then  the  voltage  across  the  plates  is  given 


VQ  =  KPL 


(19) 
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where 


2.5  Construction  of  the  Quartz  Detector  Mount 

The  perspective  view  of  the  quartz  detector  mount  which  Is 
an  assembly  of  the  crystal  holder,  rotating  mechanism  and  the  built- 
in  pre-amplifier,  is  shown  in  Figure  8.  A  metallic  cylinder  encloses 
a  plastic  cyrstal  holder  with  electrodes  and  the  preamplifier.  The 
preamplifier  is  mounted  on  the  rear  of  the  cylinder.  The  metal 
cylinder  is  supported  at  the  ends  by  two  flanges  mounted  on  a  common 
base  such  that  it  can  be  rotated  about  its  axis.  The  hcles  on  the 
front  and  rear  end  of  the  cylinder  permit  the  laser  beam  to  travel 
through  the  crystal  and  out  without  any  obstruction. 

A  cut-away  view  of  the  crystal  holder  and  electrode 
assembly  is  shown  in  Figure  9.  This  part  of  the  assembly  is  made 
with  insulator  material  to  minimize  any  disturbance  of  the  potential 
field  configuration  set  up  by  the  d.c.  Dolarization.  The  two  flanges 
with  circular  holes  at  the  center  keep  the  cylindrical  crystal  in 
position  and  do  not  interfere  with  the  path  of  the  laser  beam. 

The  primary  requirements  of  the  preamplifier  are  that 


it  should  have  a  high  input  impedance  and  low  noise  figure.  A 
balanced  cathode  follower  using  the  subminiature  tube  CK  6112  is 
built  inside  the  mount.  The  filament  is  operated  from  a  d.c. 
source  and  the  tube  is  operated  under  heavy  space  charge  conditions. 


2.6  Experimental  Arrangements 


The  general  set  up  of  the  experiment  is  shown  in  Figure  10. 


The  experiment  is  conducted  with  a  pulsed  ruby  laser  than  can 
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deliver  approximately  one  joule  of  energy  output  with  a  pulse 
width  of  approximately  400  microseconds.  The  output  from  laser 
is  passed  through  the  quartz  detector.  The  quartz  detector  mount 
has  facility  for  vertical,  transverse  and  rotational  alighments. 
The  beam  emerging  from  the  detector  strikes  a  white  background. 

The  scattered  light  is  picked  up  by  a  photomultiplier  and  is  fed 
to  one  of  the  inputs  of  the  dual  beam  oscilloscope  (Tektronix 
Model  555) .  The  output  of  the  quartz  detector  is  fed  to  the 
second  input  of  the  oscilloscope.  Both  traces  of  the  oscilloscope 
are  synchronized  with  the  trigger  voltage  that  fires  the  laser. 

The  oscilloscope  is  set  for  single  sweep  operation. 

2.7  Observation  of  D.C.  Polarization 

With  the  set  up  shown  in  Figure  10,  the  quartz  output 
wac  measured.  The  output  is  of  the  order  of  40  microvolts.  This 
is  for  the  angular  position  of  the  mount  for  which  the  output  is 
positive  maximum.  Figure  11a  is  the  oscilloscope  picture  showing 
this  output.  The  upper  trace  is  the  photomultiplier  output 
representing  the  laser  output.  The  photomultiplier  tube  circuit 
nas  a  long  time  constant  and  hence  the  spikes  in  the  ruby  output 
are  not  seen.  Only  the  envelope  of  the  laser  output  is  present. 
The  lower  trace  represents  the  quartz  detector  output  at  the 
position  of  maximum  positive  output.  Here  also  one  does  not 
see  the  spikes  though  they  are  present,  due  to  limited  bandwidth 
of  the  prea  ,f,ier.  Figure  lib  represents  the  situation  when  the 
quartz  crystal  is  replaced  by  a  glass  rod.  The  angular-  position 
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of  the  quartz  mount  Is  maintained  the  same  as  in  the  earlier  case. 
The  output  of  the  quartz  was  due  to  its  crystalline  character. 
However>  glass  is  amorphous  and  should  yield  no  output.  Thus,  a 
comparison  of  the  two  output  proves  that  the  output  voltage  i3 
due  to  the  crystalline  nature  of  the  medium.  In  Figure  lib  one 
observes  a  very  small  output  from  the  glass  rod  which  is  found 
to  be  independent  of  the  angular  position.  This  small  output 
might  be  due  to  the  residual  strain  in  the  glass  rod. 

2 . 8  Angular  Dependence  of  D.C.  Polarization 

It  was  proved  in  Section  2  that  rotating  the  quartz 
detector  about  its  axis  while  keeping  the  electric  field  orien¬ 
tation  of  the  laser  beam  fixed  in  3pace,  the  d.c.  polarization 
should  rotate  at  double  the  angular  rate.  Thus,  if  the  quartz 
detector  is  rotated  through  90°  in  either  clockwise  or  anti¬ 
clockwise  direction  the  d.c.  polarization  vector  should  rotate 
through  an  angle  of  l80°.  This  means  that  one  should  expect  a 
reversal  in  the  d.c.  output  voltage  of  the  quartz  detector  as  it 
is  rotated  u. .rough  90°.  This  is  proved  by  the  pictures  shown  in 
Figure  12.-  Figure  12a  represents  an  angular  position  for  which 
the  output  is  a  negative  maximum.  Figure  12b  represents  the 
angular  position  90°  away  from  that  corresponding  to  Figure  12a. 

It  is  seen  that  for  the  same  laser  power,  depicted  by  the  upper 
traces,  the  magnitude  of  the  quartz  output  remains  the  same  while 
the  direction  alone  reverses. 

2.9  Determination  of  Second  Order  Nonlinear  Coefficient 


It  is  possible  to  estimate  the  value  of  the  second  order 
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nonlinear  coefficient  a  in  the  matrix  of  Eq.  (6)  by  measuring  the 

d.c.  output  voltage  of  the  quartz  detector.  The  procedure  is  to 

first  determine  the  peak  power  in  the  laser  pulse  by  measuring 

the  total  snergy  using  calorimetric  technique  and  knowing  the 

pulse  shape  given  by  the  photomultiplier  output.  Knowing  this 

value  of  PL  in  Eq.  (19)  and  the  circuit  parameters  in  Figure  4, 

a  can  be  calculated  directly.  The  value  obtained  in  this  experi- 

-8 

ment  is  1.18  x  10"  e.s.u.*  In  obtaining  this  value,  correction 

has  been  made  for  the  presence  of  optical  activity  xn  quartz. 

It  can  oe  shown10  that  the  maximum  corrected  polarization  due 
to  optical  activity  will  be  0.38  times  that  of  uncorrected  value. 

It  was  mentioned  in  Section  2.1  that  the  second  order 
nonlinear  polarization  coefficient  was  related  to  the  electro¬ 
optic  coefficients.  This  has  been  worked  out  in  detail  by 
Bass12.  Following  this  approach  it  can  be  shown10  that  a  is 
given  by 

4 

n  _  q 

a  =  r^  =  0.28  x  10  e.s.u- 

where  nQ  is  the  ordinary  index  of  refraction  and  r^  the  electro- 
optic  coefficient.  Comparing  this  with  the  experimental  value, 
it  is  seen  that  the  measured  value  is  approximately  four  times 
that  of  the  theoretical  value.  This  is  attributed  to  the  errors 


*  This  value  Is  different  from  that  quoted  in  Reference  11, 
because  there  wa^  an  error  of  a  factor  of  6  in  the  initial 
measurement  of  power.  Also  the  reduction  in  polarization 
due  to  optical  activity  was  not  accounted  for  originally. 


in  measurements  of  the  exact  value  of  d.c.  polarization  and  the 

peak  power  in  the  laser. 

2.10  Splklngs  in  the  D.C.  Polarization 

The  experimental  results  of  the  previous  sections  were 

obtained  using  a  low  power  laser.  Consequently,  the  output  of 

the  quartz  detector  was  low.  The  gain-bandwidth  criteria  of  the 

detecting  system  precluded  the  possibility  of  the  observation  of 

the  spiking  phenomena  which  occurs  at  intervals  of  the  order  of 

a  microsecond  in  the  d.c.  output.  Because  of  the  interest  in 

the  feasibility  of  using  this  method  to  measure  power  in  the 

laser  beam,  a  high  power  laser  that  is  capable  of  delivering 

100  joules  output  for  500  microseconds  was  developed.  (It  should 

be  borne  in  mind  that  a  Q-switched  laser  will  not  serve  the  present 

purpose.  This  is  due  to  the  fact  that  the  Q-switched  laser 

delivers  much  higher  intensity  for  a  very  short  duration  and 

thus  the  gain-bandwidth  requirements  of  the  detecting  system 

has  the  same  restrictions.)  Now  that  the  d.c.  output  was  large, 

a  new  amplifier  using  RCA  8056  nuvistors  was  built  which  had 

higher  bandwidth  and  gain.  The  amplifier  circuit  diagram  is 

shown  in  Figure  13 .  The  8056  nuvistors  were  chosen  because  of 

their  small  size  and  electrical  properties.  They  are  low  noise, 

7 

high  frequency,  tubes  capable  of  operating  with  10  ohm  grid 
resistors.  To  achieve  gain  in  the  input  stage  and  still  avoid 
the  Miller  effect  capacitance,  a  cascode  configuration  is  used. 

The  observed  spikings  in  the  d.c.  polarization  are 
shown  in  Figure  14.  Both  in  Figures  l4a  and  14b,  the  upper 
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trace  represents  the  d.c.  output  and  the  lower  trace  the  laser 
output.  It  can  be  observed  that  the  splkings  in  the  d.c, 
polarizations  are  rounded  off  due  to  bandwidth  limitation  in 
the  detecting  system.  Also,  according  to  our  preliminary  results 
we  do  not  see  any  systematic  correlation  in  amplitudes  between 
the  fundamental  and  d.c.  splkings. 

2.11  Summary 

The  d.c.  polarization  that  is  developed  in  a  nonlinear 
dielectric  medium  of  quartz,  when  a  high  intensity  laser  beam 
propagates  through  it  is  investigated.  A  linear  relationship 
between  laser  intensity  and  d.c.  polarization  is  derived  for 
the  case  of  z-axis  propagation.  A  convenient  method  to  detect 
the  d.c.  polarization  is  suggested.  The  interaction  between 
the  propagating  wav°  and  the  detecting  circuitry  is  analyzed 
and  an  equivalent  circuit  model  developed.  It  is  proved  that 
the  phenomenon  of  d.c.  polarization  cannot  transfer  any  d.c. 
power  to  circuitry  from  the  laser  beam.  Details  of  the  experi¬ 
mental  arrangement  and  the  construction  of  the  quartz  detector 

mount  are  given.  The  second-order  nonlinear  polarization  co- 

-8 

efficient  a  of  quartz  is  estimated  to  be  1.18  x  10“°  e.s.u. 
which  is  approximately  four  times  th^  theoretical  value  predicted 
The  splkings  that  are  present  in  the  laser  can  also  be  observed 
In  the  d.c.  polarization.  However,  the  preliminary  results  do 
not  indicate  any  amplitude  correlation  between  them. 
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3.  OPTICAL  SUBHA RMCNIC  GENERATION 


3.1  Introduction 

At  the  present  time  there  is  a  serious  gap  in  signal 
sources  in  the  frequency  range,  say  300  to  5>0C0  kmc.  In 
addition  to  the  academic  interest  in  filling  the  gap  in  the 
electromagnetic  spectrum,  the  availability  of  signal  sources 
would  make  possible  some  interesting  new  spectroscopic  studies 
of  materials.  Now  that  there  is  available  for  us  very  high 
intensity  sources  at  infra-red  frequencies,  it  strongly  suggests 
the  possibility  of  achieving  this  goal  by  generating  subharmonics. 

The  above  long-range  objective  is  the  motivation  for 
the  present  study,  the  purpose  of  which  is  to  investigate  the 
feasibility  of  subharmonic  generation  at  optical  and  Infrared 
frequencies  by  interaction  of  a  high  intensity  laser  beam  with 
a  nonlinear  dielectric  medium.  ,e  principle  that  is  under 
study  is  the  parametric  subharmonic  generation. 

Interaction  of  laser  beams  with  nonlinear  dielectric 
media  has  successfully  been  employed  for  the  generation  of 
harmonics.  We  propose  here  that  a  similar  interaction  phenomenon 
will,  under  favorable  conditions,  generate  subharmonics  of 
different  orders.  Specifically,  using  a  crystal  that  lacks 
inversion  symmetry,  a  subharmonic  which  is  half  the  frequency 
of  the  laser  can  be  generated.  The  basic  principle  of  the 
phenomenon  can  best  be  explained  by  making  use  of  the  treatment 
of  parametric  oscillations  in  circuits  using  varactor  uiodes 
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(or  saturable  core  inductors).  The  treatment  presented  here  is 
of  a  general  nature,  instead  of  a  parametric  degenerate  amplifier, 
and  is  intended  to  illustrate  the  physical  nature  of  the  phe¬ 
nomenon.  In  what  follows,  the  analogy  between  an  atomic,  model 
and  circuit  theory  model  for  the  existence  of  parametric  sub¬ 
harmonic  oscillation  is  presented. 

3.2  Theoretical  Discussion 

3.2.1  Circuit  Theory  Model^ 

Consider  the  circuit  shown  in  Figure  15.  A  varactor 
(parametric  diode)  whose  capacitance  is  varied  at  an  angular 
frequency  of  2o>  is  connected  across  a  parallel  resonant  circuit. 
The  capacitance  of  the  varactor  is  given  by 

C(t)  *  C  (1  +  ~  cos  2o>t)  (21) 

°  °o 

The  resonant  frequency  of  the  circuit  with  AC  =  0  is  given  by 


“  =  rLeon/2  <22> 

If  Qe  is  the  charge  on  the  capacity  ice,  we  can  write  the  differen¬ 
tial  equation  which  describes  the  circuit  as 

0  dQ 

- i  +  C”  "hI  +  ®  U  “  ^  cos  2o>t  +  ...)  Q  =  0  (23) 

dt  o  o 

“Get/2C0 

Letting  Qg  =  e  U(t)  and  neglecting  the  higher  order  term 

inside  the  parenthesis  in  Eq.  (23),  we  obtain  the  canonical 
Mathieu  equation  which  is  of  the  following  form: 
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where 


^-2  +  (a  -  2q  cos  2a>t)  U=0 
dtT 


a 


=  1 


°e2  . 

4o)2C 

o 
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(24) 


AC 

■4c. 


(26) 


Equation  (24)  has  two  kinds  of  solutions:  stable  and  unstable. 
These  are  described  in  detail  by  Chang  .  For  the  present  dis¬ 
cussion,  it  is  enough  to  mention  that  under  unstable  conditions, 
the  circuit  sees  a  neg\tive  conductance  given  by 


Gn=-^sin20  (27) 

where  0  is  the  phase  angle  of  the  subharmonic  oscillation.  Sub¬ 
harmonic  oscillation  occurs  when  G  =  G  . 

n  e 

3.2.2  Interaction  of  Electromagnetic  Radiation  with 
Nonlinear  Dielectric  Medium 

The  phenomenon  of  the  interaction  of  radiation  with 

th^  nonlinear  dielectric  medium  will  be  illustrated  by  considering 

a  one-dimensional  enharmonic  oscillator  that  is  subjected  to  the 

forcing  field  caused  by  the  radiation.  In  the  absence  of 

5 

radiation,  the  Hamiltonian  of  the  system  is  of  the  form 

H  -  |e  +  ™02x2  -  3  *  x3  -  (28) 


where  m  is  the  mass,  p  the  momentum  of  the  particle,  X  and  T| 
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constants  of  the  system,  and  x  is  the  displacement  of  the  particle 
from  its  equilibrium  position.  In  the  electric  dipole  approxima¬ 
tion,  the  interation  of  the  system  with  radiation  is 


Hlnt  «  exEQ  cos  (2o)t  +  0)  (29) 

where  e  is  the  electronic  charge  and  Eq  cos  (2cnt  +  0)  is  the 
incident  monochromatic  radiation.  Prom  Eqs.  (28)  and  (29)  the 
equation  of  motion,  including  a  dissipative  term,  for  the 
system  can  be  written,  which  will  assume  a  general  form  as 
below 

2 

+  a  4-  +  P  x  +  y  x2  +  P  x3  =  B  cos  2a>t  (30) 

dt* 
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Following  the  procedure  of  Hayashi  ,  assuming  a  solution  of 
the  form 


x(t) 

we  can  obtain 


=  xo  +  kl  sin  ^  +  k2  C0S  ^  +  k3  COS  2ait 
from  Eq.  (30)  the  following  Mathieu  equation 


(31) 


cos(  vcct 


(32) 


V -1 

Hayashi1^  discusses  the  stability  conditions  of  Eq.  (32).  Under 

unstable  conditions,  the  system  will  oscillate  at  half  the 

radiation  frequency,  thus  generating  subharmonics.  We  can  also 

recognize  this  by  comparing  Eqs.  (24)  and  (32). 

3.2.3  Discussion  on  the  Analogy  between  Circuit  and 
Field  Theory  Concepts 

The  above  discussion  mathematically  portrays  the 
similarity  oetween  subharmonic  generation  in  circuitry  and 
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atomic  systems.  This  can  also  be  explained  qualitatively  by 

considering  the  physical  mechanisms.  In  the  circuitry  model. 

the  varactor  capacitance  is  varied  at  the  pump  frequency.  Tnis 

can  be  visualized  as  the  constant  of  the  dielectric  in  the 

capacitor  being  varied  at  the  pump  frequency.  In  the  case  of 

interaction  of  radiation  with  the  dielectric,  the  dielectric 

constant  is  varied  at  the  radiation  frequency  due  to  the  second 

order  nonlinear  term.  That  is, 

e  =  e  +  a  E  cos  ait 
o  o 

where  e  is  the  total  dielectric  constant,  e  the  dielectric 

o 

constant  due  to  linear  term  between  polarization  and  electric 
field  intensity,  a  constant  of  proportionality  and  Eq  cos  cut 
the  radiation  field  intensity. 

Although  we  have  attempted  to  show  the  similarity 
between  circuit  model  arid  field  theory  model,  we  should  bear 
in  mind  the  limitations  of  extending  the  concept  of  circuit 
theory  into  the  optical  region.  The  availability  of  high  power 
monochromatic  sources  at  optical  and  Infrared  frequencies  was 
a  major  problem  till  now.  This  has  been  alleviated  with  the 
advent  of  lasers.  The  second  Impediment  is  the  proper  phase¬ 
matching  between  the  fundamental  and  subharmonic  so  as  to 
have  an  efficient  transfer  of  power.  The  problems  involved 

here  are  identical  to  those  of  the  generation  of  harmonics 

2  5 

and  have  been  discussed  elsewhere  in  detail  '  .  It  can, 
however,  be  pointed  out  in  our  discussion  cf  the  analogy, 
that  the  phase-matching  condition  is  easily  solved  in  the 
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circuitry  by  the  presence  of  the  idler.  In  the  field-interaction 
case  the  extraordinary  wave  of  one  frequency  is  matched  to  the 
ordinary  wave  of  the  other. 

3.3  Experiment 

The  experimental  arrangement  is  shown  in  Figure  16. 

The  output  of  the  laser  J s  a  collimated  beam  with  a  maximum 
energy  of  approximately  50  Joules.  In  addition  to  the  laser 
light  the  output  contains  visible  and  infrared  radiation  from 
the  flashtubes. 

Since  the  signal  to  be  detected  has  a  wavelength  of 
1.39 \i,  the  radiation  from  the  flashtubes  in  this  range  must  be 
filtered  out  prior  to  the  anisotropic  cyrstal.  A  four- inch 
cell  filled  with  glycerin  is  used  for  this  purpose.  A  trans¬ 
mission  curve  for  1.1  cm.  of  glycerin  (in  a  plexiglass  cell) 
is  shown  in  Figure  17.  As  can  be  seen  from  this  curve,  the 
attenuation  at  69^3  A°  is  negligible  while  the  attenuation  at 
1.39|i  is  very  large,  approximately  93#. 

The  attenuation  increases  exponentially  with  length. 
Hence,  the  four-inch  cell  attenuates  radiation  with  a  wave¬ 
length  greater  than  l,37p-  by  at  least  60  db.  The  maximum 
transmission  for  wavelengtns  greater  than  1.2p,  occurs  at  1.32;x, 
and  even  this  is  attenuated  by  more  than  30  db.  Thus,  the 
primary  signal  leaving  the  glycerin  cell  is  the  laser  beam  at 
69^3  A° . 

If  the  crystal  Is  such  that  it  will  support  the 
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subharmonic,  then  the  output  will  consist  of  the  laser  beam 
(very  intense)  plus  the  subharmonic  of  the  laser  beam.  In  order 
to  detect  the  subharmonic  signal,  the  laser  beam  at  69^3  and 
all  other  radiation  outside  of  a  small  band  at  1 . 39m-  must  be 
suppressed.  This  is  accomplished  by  a  silicon  filter. 

The  band  gap  energy  in  pure  silicon  is  1.1  e.v.; 
thus  silicon  has  a  very  high  absorption  coefficient  for  radiation 
with  a  wavelength  of  l.lp,  or  less.  The  absorption  curve  is  shown 
in  Figure  18.  From  measurements  made  on  a  10-mil  wafer  of 
silicon,  the  attenuation  of  the  laser  beam  in  the  3  mm  thick 
disc  used  should  be  greater  than  180  db. 

The  use  of  silicon  has  disadvantages  in  that  the  index 
of  refraction  at  1.39p.  is  3-97,  which  means  the  reflection  co¬ 
efficient  is  .33.  In  addition  to  this,  our  silicon  disc  is 
slightly  wedge-shaped  and  also  apparently  acts  as  a  scatterer. 

The  net  result  is  an  8  db  insertion  loss  at  1.39u* 

The  detection  is  accomplished  by  using  a  PbS  photo- 
conductive  cell  mounted  on  a  three-way  positioner  so  that  it 
may  be  positioned  at  the  point  of  maximum  signal.  The  signal 
waveform  i3  observed  on  a  Tektronix  355  oscilloscope  with  a 
Type  E  head. 

*.s  a  final  prevention  against  unwanted  radiation,  a 
dielectric  interference  filter  will  be  mounted  directly  in 
front  of  the  PbS  device. 
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Axes  Orientation  for  Deriving  the  Angular  Dependence 
of  d.  c.  Polarization  for  Propagation  Along  z-axis 


Simple  Parallel  Plate  Capacitor  Formed  with  the 
Nonlinear  Dielectric.  The  Polarization  P  is  the 
External  Polarization  Induced  by  the  Laser  Beam. 


Concentri 


Equi potential  Linen  Outside  the  Quartz  Medium 


SUPPORTING  FLANGES 


Perspective  Viev  of  the  Crystal  Moult 


Cut-Away  View  or  Quartz  Detector 


(a) 


/«.  \ 
\U> 


Ti g.  11  Comparison  of  Output  from  Quartz  Crystal  with  that  fran 

Glass  Rod 

(a)  Photograph  of  the  duel  beam  oscilloscope  in  which 
the  lover  trace  shows  the  quartz  detector  output 
and  the  upper  trace  the  intensity  of  the  laser 

beam. 

Sweep  rate:  5  00  microseconds/cm. 

Sensitivity  of  upper  trace:  50  raicrovolts/em. 

(b)  As  in  (a)  with  the  quartz  replaced  by  the  glass 
rod. 


(b) 


Fig.  :.2  Angular  Dependence  of  d.  e.  Polarization 

(a)  Photograph  of  tie  dual  bean  oscilloscope  in  which 
the  lower  trace  shows  the  quartz  detector  output 
and  the  upper  trace  the  intensity  of  the  laser 
beam.  Quartz  detector  Is  adjusted  for  maxitmsn 
negative  output. 

(b)  As  in  (a);  quartz  detector  rotated  by  90° 
from  that  of  position  in  (a). 


. . . 


FIG. 13  AMPLIFIER  USED  TO  MEASURE  SPIKINGS 
IN  THE  DC  POLARIZATION. 


Pig.  14  Spikings  in  the  Polarization 

(a)  Photograph  of  the  dual  beam  oscilloscope 
in  which  the  lower  trace  shows  the  d.c. 
polarization  and  the  upper  trace  the  laser 
output.  Time  scale:  5  micro3econds/cm. 

(b)  As  in  ^a)  with  the  time  scale:  2  micro¬ 
seconds/  cm. 


0 

c 

[ 

[ 


FIG. 15  BASIC  CIRCUIT  FOR  SUBHARMONIC  GENERATION. 
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